Stability of half quantum vortex in rotating superfluid ^He-A between parallel plates 
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We have found the precise stabihty region of the half quantum vortex (HQV) for superfluid '^He 
A phase confined in parallel plates with a narrow gap under rotation. Standard Ginzburg-Landau 
free energy, which is well established, is solved to locate the stability region spanned by temperature 
T and rotation speed (Q). This Q-T stability region is wide enough to check it experimentally in 
available experimental setup. The detailed order parameter structure of HQV characterized by Ai 
core is given to facilitate the physical reasons of its stability over other vortices or textures. 
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Half quantum vortex (HQV) and associated Mojorana 
zero energy mode have been widely discussed in vari- 
ous research fields in condensed matter physics, ranging 
from superconductors, superfluids, graphene and frac- 
tional Quantum Hall systems[l]. In particular, theoreti- 
cal and experimental investigations are devoted to finding 
HQV in superconductors and neutral Fermion superflu- 
ids in cold atoms. Recently strong interest on HQV is 
partly motivated by the fact that the bound state cre- 
ated in the core of HQV is characterized by the Majo- 
rana state with the zero energy exactly at the Fermi level. 
The Majorana particleQ is thought to be a candidate 
for quantum computation because it obeys non-Abelian 
statistics and its existence is protected topologically to 
avoid decoherence. These situations are ideal for quan- 
tum computation!^ if it really exists. 

So far there has been no firm experimental evidence for 
HQV in any superconductors. It is necessary for HQV 
to exist that superconductivity is described by a chiral 
p-wave pairing where the d vector is able to be free to 
rotate. It has often been argued that Sr2Ru04 may be a 
prime candidate 0, |3, 0], but strong doubt has been cast 
on this possibility of Sr2Ru04 of its triplet pairing[! 0, 
[loj . Note that the first discovered triplet superconductor 
UPta is an /-wave pairing, not chiral p-wave 11 1. 

Superfluid ^He-A phase is characterized by a chiral 
wave pairing. There is no doubt on this identification!. 
In fact, Volovik and Mineev[l^ are the first to point out 
the possibility to the realization of HQV in 1976. Since 
then, there have been several general arguments on the 
stability of a HQV in connection with ^He-A phase [13. 
[H, [l3|. However, there are no serious calculations which 
consider realistic situation in superfiuid '^Hc-A phase on 
how to stabilize it and on what boundary conditions are 
needed for it. 

Recently, Yamashita, et al[l3l have performed an ex- 
periment intended to observe HQV in superfluid '^He-A 
in parallel plate geometry. The superfluid is confined in 
a cylindrical region with the radius R = 1.5mm and the 
height 12.5/im sandwiched by parallel plates. A mag- 
netic field H — 26.7mT(||z) is applied perpendicular to 
the parallel plates under pressure P=3.05MPa. Since 



the gap 12.5/im between plates is narrow compared to 
the dipole coherence length ^ 10/J,m, the / vector is 
always perpendicular to the plates. Also the d vector is 
confined within the plane because the dipole magnetic 
field Hd ^ 2.0mT. They investigate to seek out various 
parameter spaces, such as temperature T, or the rota- 
tion speed f2 up to r2 = 6.28rad/s by using the rotating 
cryostat in ISSP, Univ. Tokyo, capable for the maxi- 
mum rotation speed ~12rad/s, but there is no evidence 
for HQV 17| . Here we are going to give an answer why 
it is so and to examine the stability region of a HQV 
which competes with the ordinary singular vortex with 
the integer winding number and propose a concrete ex- 
perimental setup which is feasible to perform in the light 
of the present experimental situation[17]. 

We start by examining the possible order parame- 
ter (OP) forms allowed under the above experimental 
conditions. The OP of the superfluid '^He is given by 
A = «Sp,i(v4^,icrpPi)(Tj, in general where A^^i isa3x3ma- 
trix (/i, i — x,y, z). Among the two known bulk phases as 
ABM (A) and BW (B) phases, we focus on the A phase 
in this paper, which is written as ^ = Aorfp(r? + irh) 
conveniently expressed in terms of the d vector and the 
I vector [l — n X rh, n and rn are unit vectors forming a 
triad). The former (latter) vector characterizes the spin 
(orbital) part of the OP. Since the / vector is locked per- 
pendicular to the plates along the z axis and the applied 
fleld {H II z) perpendicular to the plates confines the d 
vector within the plane (a;, y) because it is strong enough 
compared with the dipole field as mentioned above. 

The HQV form originally proposed by Salomaa and 
Volovik [14], which was followed by others]! d, 0] , can be 
expressed as 



irn)i 



(1) 



where d = x cos sin | {9 is the angle from the x axis) . 
Since the d vector is assumed to be real here, we call it 
R-HQV. When winding around the vortex core, the OP 
exhibits simultaneous change of sign of the d vector and 
shift of the phase 9. Namely, (6*, d) => [9 + tt, -d). The 
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TT phase windings of the orbital and the spin parts add 
up, resulting in 27r phase winding in total. Alternatively, 
the wave function of this R-HQV form is cast in a form 



V-^ Ao(r)(e^^| Tt) + Ui))(p. + *Pj,). 



(2) 



It is clear that in the R-HQV the pairs phase- wind 
by 27r while the J.| pairs do not. It will turn out shortly 
that this somewhat restrictive R-HQV form is not a full 
solution of our Ginzburg-Landau (GL) free energy func- 
tional under rotation. Thus we have to seek more general 
HQV solution to be competitive with the vortex free state 
stable at rest (A phase texture; AT) and the singular vor- 
tex with integer winding (SV). We generalize the OP to 
find the stable HQV by noticing that the orbital part is 
doubly degenerate p ± ipy — p± in addition to doubly 
degenerate spin space. The most general wave function 
is spanned by four basis functions, namely, 

^^{A++ir,e)\^^)+A^+ir,e)\U))p+ 

+ {A+_{r,9)\ tT>+A__(r,0)| U))p_ (3) 

where each component A±±{r,9) = A±±(r)e*''™±± can 
have its own winding number 'w±± in the polar coordi- 
nates. Under axis-symmetry w±+ = w±- — 2n must be 
satisfied with n > being integer [l8|. The winding num- 
ber combination {w^+jW i-,w^ — (1,0,3,2) is 

straightforwardly generalized from the above R-HQV 
form (2), which we are trying to stabilize. The (0,0,2,2) 
phase is the A phase texture (AT) and (1,1,3,3) is the 
ordinary singular vortex (SV). The AT is always stable 
at rest, and HQV and SV compete each other under ro- 
tation. Other several phases with different winding num- 
ber combinations, such as (0,-1,2,1), (-1,-2,1,0) or (-1,- 
1,1,1) are all irrelevant; namely they are never stabilized. 
Note that the (0,-1,2,1) phase is stable next to the lowest 
AT (0,0,2,2) at rest. Under counter clock-wise rotation 
the chiral p+ ^ Px + iPy is favored over p- = p^ — ipy 
Thus the p+(p_) component constitutes the major (mi- 
nor) one. 

The Ginzburg-Landau (GL) free energy functional is 
well establishedilj, [H, M, 0, EH and given by a stan- 
dard form 
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FIG. 1: (Color online) Free energy comparison for AT, HQV, 
SV and (0,-1,2,1) states as a function of Q. for R = 10/im and 
t = T/Tc = 0.95. 
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where di — Vi 

aoil-T/T, + f,AT/T,) (fi = ±, ao = ^), 
7C(3)iV(0)(nui.)V240(7rfcBTc)2. is the coupling con- 
stant of the dipole interaction, which is gd <IC an[l3|. 
As mentioned above, we assume a two-dimensional sys- 
tem for the OP spatial variation; fi,i = x,y or ±. 
The magnetic field acts not only to pin the d vec- 
tor within the plane, but also to shift the transition 
temperature by At = AT/Tc = {T^i - Tct)/2Tc. 
The fourth order GL coefficients are given by /3i — 
-(1 + 0.1<5)/3o, /32 = (2 + 0.2<5)/3o, /Sa = (2 - 0.05<5)/3o, 
/34 = -(2 - 0.0555)/3o, and /Sg = -(2 + 0.15)(iQ where 
/3o = 7C(3)iV(0)/120(7rfcBTc)^[il]. The strong coupling 
correction (5 > due to spin fluctuations serves sta- 
bilizing the A phase over the B phase in the (P, T) 
phase diagram [2^. In the following we use the GL 
parameters fi^ tabulated lol. 2oL 21 1 appropriate for the 
experiment at P=3.05MPa. 

We find the free energy minima under the rigid bound- 
ary condition = for r > i? (i? is the radius of the 
system). A fundamental difficulty associated with the 
numerical computations lies in the fact that the coherent 
length ^ = lOnm is extremely small compared with the 
system size R where we have to take care of these two 
length scales in the equal footing in order to accurately 
evaluate the relative stability among three textures; AT, 
HQV and SV. This is a reason why this kind of serious 
energy comparison has not been done before. We care- 
fully calibrate the accuracy of our numerical computation 
to allow the detailed comparison. 

We first consider the weak field case where the transi- 
tion temperature splitting At ~ 0. As shown in Fig. 1 we 
compare three phases with additional other phase men- 
tioned above. It is seen that at rest and lower rotation 
region AT is stable and eventually upon increasing Jl, 
SV takes over at f2c2- Although the HQV is stabler than 
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FIG. 2: (Color online) Order parameter and d vector profiles 
for AT (a), HQV (b) and SV (c) for R = 1.0/im and t = 0.97. 
Left and center columns show the cross-section of OP along 
the radial direction r. Right column shows d vector patterns. 
In (b) it winds by n around the center r = 0. 



SV at rest situated almost at the half way between AT 
and SV because the phase winding occurs only for the 
pairs. Under rotation the energy gain due to the angular 
momentum is less than that in SV because of the above 
reason. Thus HQV is never stabilized under weak field 
region. We also plot the (0,-1,2,1) state which is second 
lowest at rest and becomes irrelevant under rotation. 

Note that the previous form (2) of R-HQV does not im- 
prove this situation, rather becomes worse in its stability. 
The R-HQV (2) indicates that the vortex core singularity 
occurs for both and [[ pairs even though the latter 
does not have phase-winding, leading to the additional 
loss of the condensation energy. The strong coupling ef- 
fect acts to destabilize both R-HQV and SV relative to 
AT, thus R-HQV is never stabilized (see below). 

In Fig. 2 we illustrate the results of the OP profiles for 
AT (a), HQV (b) and SV (c) and their d vector textures. 
The OP's in AT are uniform in the central region around 
r — Q, decreasing towards zero at the boundary r ^ R 
whose characteristic length is ^. Thus AT is basically 
the A phase in the bulk. Near the boundary the induced 
components and appear peripherally. 

In HQV (Fig. 2(b)) one of the two majority compo- 
nents with = 1 exhibits a phase singularity at 

r = 0, the other component A y with w y = being 

depressed slightly there. A^ with = 3 and A 

with w = 2 are also induced at r = and r — R. 

Therefore this HQV profile shows that only the jj pairs 
appears at around r = 0, implying the Ai core state. 
This tends to stabilize this HQV further compared with 
R-HQV given by Eq.(2) because (A) the condensation 
energy loss is less, (B) the fourth order GL energy con- 
cerning the interaction term between |f and || pairs can 
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FIG. 3: (Color online) Stability region of HQV sandwiched 
between fici and Qc2 as a function of R {t — 0.97 and 
At = 0.05). f2ci=0.05 rad/sec is the extrapolated value 
for i?= 1.5mm. Inset shows the free energy comparison for 
R — 20^m, displaying the successive transitions from AT to 
HQV at Q.C1 and from HQV to SV at r2c2. 



be expressed as — 4(5/3o|d+p|(i_p. This particular term 
due to the strong coupling acts to earn the extra gain for 
this HQV. However, AT is simultaneously stabilized by 
this term, thus HQV never wins in weak fields. Note in 
passing that the IT ^ind J,| pairs are completely indepen- 
dent when (5 = because the weak coupling GL form is 
derived under the assumption that the spin space is rota- 
tionally invariant. It is seen from Fig. 2 that the d vector 
rotates by n when going around the origin in HQV (b) 
while in the others (a) and (c) it is uniform. 

Finally SV in Fig. 2(c) exhibits the phase singularities 

for both major components j4++ with w++ = 1 and A ^ 

with w ^ ~ 1 and the induced components A^ with 

W-i — 3 and A with w = 3 appear at the places 

where the OP spatially varies. Thus this SV is quite 
advantageous under rotation because they can absorb ef- 
ficiently the rotational kinetic energy. 

Having found that HQV is not stable in weak field re- 
gion {H ~ Hd =10mT) both at rest and under rotation, 
we resort to higher field region; an order of a few kG 
where At ^ or Tc| ^ . This extension indeed stabi- 
lizes the HQV as shown in the inset of Fig. 3 where we 
compare the three states as a function of Q. It is seen 
that as increasing il, AT changes into HQV at fici and 
then HQV to SV at ^lc2- The relative stability region 
flci/Qc2 ^ 0.8 which is wide enough to check experimen- 
tally. 

The reason for the HQV stabilization is physically ex- 
plained as follows. (A) By introducing At which increases 
(decreases) the OP amplitude of jj pair (ft pair), the 
kinetic energy loss due to the phase winding of A^^ 
with = 1, which remains unscreened and spreads 

out whole system, becomes less compared to AT at rest, 
meaning that the HQV energy approaches towards the 
AT energy in Fig. 1 at = as seen from inset of Fig. 
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FIG. 4: (Color online) Stability region of HQV in fl versus 
T/Tc {R = and At = 0.05). AiT denotes the Ai phase 

texture where only jj, pairs exist. 



3. Under rotation the HQV energy decreases by absorb- 
ing the rotation kinetic energy and eventuaUy becomes 
lower at fid , which is smaller than flc2 , stabilizing HQV 
over SV. 

The main panel in Fig. 3 shows f2ci and Qc2 as a 
function of the system size R. It is seen that the rela- 
tive stability region flci/flc2 ~ 0.8 stays at a constant 
against R, keeping 20% region above the critical rotation 
speed flci at which single HQV is created in the sys- 
tem. The extrapolated Qci to i?= 1.5mm, by which Ya- 
mashita, et aUl7| have performed experiments, amounts 
to r^ci ~ 0.05rad/sec. The rotation speed stability of 
the rotation cryostat at ISSP, Univ. Tokyo is accurate 
enough to perform it. We also notice that by changing 
the radius R of the system one can control the Qci value 
at will. For example, in i? = 300/im, fid ~ Irad/sec 
which is convenient speed to run their rotating cryostat. 

In Fig. 4 we depict the temperature dependence of 
the HQV stability region. We see that the stability re- 
gion for HQV is widen divergently when approaching the 
lower critical temperature Tc| from below where the dis- 
parity of the OP amplitudes between ft pair and || pair 
increases. Note that the actual lower transition temper- 
ature Tc| is shifted slightly upward because the spatially 

varying ft pair OP ^ induces j4_|_+(r). In other 

words, the Ai phase for > T > Tc-\ becomes narrower. 
Thus one needs not only careful temperature control, an 
order of O.OIK which is feasible enough, but also theo- 
retical backup to estimate this shift in order to precisely 
locate the HQV stability region under actual experimen- 
tal setup. 

Since the HQV has the odd winding number for the H 
pairs, the Majorana quasi-particle with zero energy ex- 
actly at the Fermi level, which is localized in the vortex 
core, is ensured by both the index theorem based on topo- 
logical argument [2J], or directly solving the Bogoliubov- 
de Gennes equation [2^. These arguments are based on 
the assumption that the ft pair and J,j pair are com- 
pletely decoupled. Here the situation is more subtle. The 
tt pair and || pair are interacting through the fourth or- 
der GL terms, which is mentioned above. These terms 



comes from the strong coupling effect due to ferromag- 
netic spin fluctuations [2^ . which ultimately help stabi- 
lizing the present HQV. Therefore, it is not obvious com- 
pletely that the present HQV can accommodate the Mo- 
jorana fermion with the exactly zero energy. This issue 
belongs to a future problem. 

We also remark on the experimental point that the 
identification of the HQV is not an easy task. The HQV 
and SV are indistinguishable by the usual NMR method 
which utilizes the satellite position in the spectrum 17 1 
because d ± I are always kept for both vortices, giving 
rise to the identical NMR spectra. We suggest small 
tilting of the field direction from H \\ z may yield the 
different NMR signatures. This point deserves further 
elaboration. 

Finally it should be pointed out that our previous the- 
ory for the parallel geometry of the superfluid ^He 25| 
differs in the field orientation H J. z there. The singu- 
lar vortex with odd integer winding number was found 
in this spinless chiral superfluid. This also gives rise to 
the Majorana zero energy mode. Thus the field orienta- 
tions yield different vortices, but those accommodate the 
Majorana particle localized at each vortex core. 

In conclusion, we have found the stability region of half 
quantum vortex in T-fl plane of superfluid '^He A phase 
confined in parallel plates and given physical reasons why 
it is stabler than A phase texture or ordinary singular 
vortex. We propose a concrete experimental setup, which 
is feasible by using the rotating cryostat such as in ISSP, 
Univ. Tokyo. 
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